State  Variable  Analysis  of  a Boiler  System 


SECURITY  CLASSIFICATION  OF  THIS  RACK  (Whmt  Dmm  Knffd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


4.  TITLE  fan*  SuOMII.j 


•erioo  covered 


State  Variable  Analysis  of  a Boiler  System.  / jMaster's  ^hesijj 


• ■ PERFORMING  ORG.  REPORT  NUMBER 


(ji 


7.  AUTHORS 

Chusakdi / Senanikrom 

J 

k 

S.  PERFORMING  ORGANIZATION  NAME  AND  ADORESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


II.  CONTROLLING  OFFICE  NAME  ANO  AOOREM 

Naval  Postgraduate  School 
Monterey,  California  93940 


(M 


lU'W.I  I-VTIE-TW 


Mar  Eli  M78 


ii.  numrer  of  faoes 

106 


14.  MONITQRINO  AGENCY  NAME  4 AOOREM (II  MIJ»^m(  ltmm  C&Mmtttnt  BfHem)  l».  SECURITY  CLASS.  (»l  thla  rmpmri) 

UNCLASSIFIED 

Naval  Postgraduate  School  y ' °^/  ^ h*  * ; 

Monterey,  California  93940  J | 


ASSIFICATI  ON/ DOWNGRADING 


IS.  OISTRI GUTION  STATEMENT  (•!  Ml*  *M«IJ 

Approved  for  public  release;  distribution  unlimited. 


17.  OISTRIRUTION  STATEMENT  (ml  (A.  AMnel  an tmrm4  in  Btmmk  JO,  II  mUrnrmt  trmm  t) 


IS.  KEY  WOROS  rCMIlim  an  kwh  mIBm  II  nmmmm 


• «4  ISnllly  If  klmm*  I 


20.  ABSTRACT  (CmmMtmtm  an  tmmmrmm  «M*  

The  state  variable  formulation  of  a Foster  Wheller  ESD-III  boiler  is 
developed  from  fundamental  principles.  The  response  of  the  model  for 
various  input  signals  is  determined  using  CSMP-III,  the  IBM  simulation  


50  i jan  7i  U73  EDITION  OF  I NOV  SS  IS  OBSOLETE 
(Page  1)  s/r  oioi-oi4*sso 1 1 


SECURITY  CLASSIFICATION  OF  THIS 


tu  cumrv  CL*»»iyiC*Tiow  or  tmh  n.<«  tnift 


^language.  The  sensitivity  of  the  model  to  various  coefficient  values  is  noted 
as  are  the  characteristics  of  various  system  states  for  small  perturbation 
values. 


Approved  for  public  release;  distribution  unlimited 


STATE  VARIABLE  ANALYSIS  OF  A BOILER  SYSTEM 


by 


Chusakdi  Senanikrom 
Lieutenant,  Royal  Tnai  Navy 
B.S.,  Royal  Thai  Naval  Academy,  iy67 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1978 


Author: 


C&MAaJidU  . 


Approved 


TEalrman,  UeparTment“of“7!ecHInicaT“Sngin€ering 


3 


ABSTRACT 


The  state  variable  formulation  of  a Foster  Wheeler 
ESD-lboiler  is  developed  from  fundamental  principles. 
The  response  of  the  model  for  various  input  signals  is 
determined  using  CSMP  -I, the  IBM  simulation  language. 
The  sensitivity  of  the  model  to  various  coefficient 
values  is  noted  as  are  the  characteristics  of  various 
system  states  for  small  perturbation  values. 
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NOMENCLATURE 


A*  .A*,  A#  = 


C 


C9  = 
Ci  = 


f»)  f*>  f » ~ 


Ka*  = 

^9S»  = 


k3»,k 


(M 


superheater,  riser  and  downcomer  cross-sectional 

2 

areas  respectively  (ft  ) 

evaporation  level  proportionality  constant 
(ft-sec/lb) 

heat  capacitance  for  superheater,  riser  and  downcomer 
tubes  respectively  (BTU/lb-8) 

average  heat  capacitance  of  combustion  gases  times 
air- fuel  ratio  (BTU/lb-R) 
neat  capacitance  for  feedwater (BTU/lb-R) 
average  heat  capacitance  of  superheated  steam 
(BTU/lb-R) 

frictior  coefficients  for  superheater,  riser  and 

2 5 

downcomer  tubes,  respectively  (sec  /ft  ) 

2 

acceleration  due  to  gravity  (ft/sec  ) 
enthalpy  of  saturated  vapour  corresponding  to  P 

B 

(BlU/lb) 

enthalpy  of  saturated  liquid  corresponding  to  P 

B 

(BTU/lb) 

enthalpy  of  liquid  in  mud-drum  (BTU/lb) 

enthalpy  of  drum  and  downcomer  liquid  (BTU/lb) 

enthalpy  of  evaporation  corresponding  to  ? (BTU/lb) 

B 

air  cooling  coefficient  at  superheater  bank  (BTU/sec) 
air  cooling  coefficient  at  riser  bank  (BTU/sec) 
heat-transf er  coefficients  from  combustion  gas  to 
superheater  tubes,  and  from  superheater  tubes 
to  steam,  respectively  (BTU/lb-*fi) 
heat-transfer  coefficients  from  combustion  gas  to 
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riser  tubes  and  from  riser  tubes  to  boiling 

• • 4 • 

liguid  respectively  (BTU/lb-H) , (BTU/lb-F  ) , (BTU/lb-R) 
Kc  = evaporation  rate  constant  of  drum  liquid  (lb/sec-R) 
Kt,K6=  constants  for  state  equations  of  saturated  steam 

o * 

( R-f  t /lb)  , (ft  | 

'-»)Lv>ld  = superheater,  riser  and  downcomer  tube  lengths 
respectively  (ft) 

M = mass  of  drum  liquid  (lb) 

Ms  = mass  of  superheater  tubes  (lb) 

M8  = mass  of  riser  tubes  (lb) 

2 

Pa  = drum  pressure  (lb/ft  ) 

2 

Ps  = superheater  outlet  pressure  (lb/ft  ) 

2 

P*  = mud-drum  pressure  (lb/ft  ) 

as = heat-input  rate  from  tube  walls  into  the  superheated 
steam  (BTU/sec) 

ft3S=  heat-input  rate  from  hot  gasses  into  superheater  tube 
walls  (3TU/sec) 

heat-input  rate  from  riser  tube  walls  into  boiling 
liquid  (BTU/sec) 

Qge=  neat-input  rate  from  hot  gasses  into  riser  tube  walls 
(BTU/sec) 

Ti  = feedwater  temperature  (R) 

tb  = saturation  temperature  corresponding  to  P (R) 

B 

T$  = superheater  outlet  temperature  (#R) 

Tw  = drum  and  downcomer  liquid  temperature  (*R) 

Ts*=  superheater  tube-wall  temperature  (°R) 

T.w=  riser  tube-wall  temperature  (*R) 

Tj*=  average  gas  temperature  at  superheater  banks  (*R) 
average  gas  temperature  at  riser  banks  (*R) 

Te=  comDustion  gas  temperature  entering  superheater  banks 
<*R) 

3 

volume  of  vapor  phase  in  drum  (ft  ) 
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= velocity  of  riser  mixture  (ft/sec) 

= velocity  of  downcomer  water  (ft/sec) 

3 

= total  drum  volume  (ft  ) 

= steam  mass-flow  rate  at  the  superheater  outlet 
(lfc/sec) 

= fuel  mass- flow  rate  (lb/sec) 

= feedwater  mass-flow  rate  (lb/sec) 

= downcomer  mass-flow  rate  (lo/sec) 

= riser  mass- flow  rate  (lb/sec) 

.=  air  mass-flow  rate  from  blower  (lb/sec) 
c=  chemically  correct  +50X  excess  air  rate  (lb/sec) 

= mass-evaporation  rate  from  drum  liquid  surface 
(lb/sec) 

= steam  mass-flow  rate  from  drum  into  superheater 
(lfc/sec) 

= quality  of  mixture  leaving  riser 
= throttle  opening  (£) 

= drum  liquid  level  (ft) 

3 

,=  saturated  vapor  density  corresponding  to  P (lb/ft  ) 

B 

3 

= superheater  outlet  density  (lb/ft  ) 

3 

= saturated  liquid  density  corresponding  to  P (lfc/ft  ) 


= density  of  liquid  vapor  mixture  leaving  riser 
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INTRODUCTION 


Boilers  as  understood  by  marine  engineers  are  closed 
vessels  containing  water  which  by  the  application  of  heat  is 
converted  into  steam  at  any  designed  pressure.  This  steam 
is  than  used  for  the  production,  through  machinery,  of 
useful  work.  A dynamic  model  of  steam  turbine  machinery 
consists  of  a boiler  model  and  a turbine  model.  Tne 
difficult  part  is  the  boiler,  as  a load  change  causes 
variations  in  some  important  properties  such  as  toiler 
pressure,  temperature  and  drum  water  level. The  analysis  in 
this  paper  will  use  Chieris  dynamic  analysis  (1)  of  a boiler 
as  a ref erence.Chien  considered  a naval  boiler  which  for 
purposes  of  analysis  was  divided  into  four  sections  namely  a 
superheater,  a downcomer- riser  loop,  a drum  and  a gas  path. 

The  principles  of  thermodynamics, heat-transfer  and  fluid 
mechanics  were  used  to  describe  the  dynamic  behaviour 
corresponding  to  each  section  of  the  boiler  and  these  were 
derived  from  equations  of  continuity,  energy, 
heat-transfer , and  momentum.  The  eguations  involve  partial 
differentials  as  well  as  nonlinear  terms.  These  eguations 
were  reduced  to  the  ordinary  linear  eguation  form  by 
applying  small  perturbation  and  difference  equation 
techniques. Linear  equations  thus  obtained  were  reduced  to 
ten  state  variable  equations  and  solved  by  digital  computer 
techniques. 

Since  there  is  an  increasing  interest  in  toiler 
modelling,  the  objective  of  this  thesis  was  to  develop  a 
comprehensive  boiler  simulation  model  in  a form  useful  for 
modern  control  (i.e.,  multivariable  control)  analysis. 
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II.  BOILER  MODEL 


A.  BOILER  CONSIDERATIONS 


The  control  problems  of  high  pressure  boilers  have 
become  more  and  more  critical,  both  from  the  operational  and 
the  economical  points  of  view.  A dynamic  analysis  is  the 
method  to  be  used  for  a control- system  analysis.  It 
consists  of  a complete  understanding  of  the  process  to  be 
controlled  and  the  effects  of  physical  and  chemical  changes. 
The  analysis  is  not  exact  by  any  means  but  the  results 
obtained  should  be  in  good  qualitative  agreement  with  actual 
tests.  The  major  difficulty  in  boiler  analysis  is  the  fact 
that  the  whole  system  is  very  complex  and  contains  numerous 
variables. 

Chien(l)  considered  a naval  boiler  which  for  purposes  of 
analysis  was  divided  into  four  distinct  sections  namely  a 
superheater,  a downcomer-riser  loop,  a drum  and  a gas  path. 
The  detail  of  analysis  is  described  in  Part  B of  this 
section.  The  basic  equations  used  in  the  analysis  are  those 
of  continuity, energy (heat-transfer) , momentum  and  the  state 
equations.  These  equations  involve  partial  differentiation 
as  well  as  nonlinearities.  Generally, the  equations  have  the 
form 

f C x , y , 2> ) > o 

. \ 

To  eliminate  the  nonLinearities  one  uses  perturbation 
theory,  which  effectively  approximates  the  response  of  the 
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system  to  small  signal  changes  about  a chosen  operating 
condition.  Thus, the  equation  is  perturbed  about  its  steady 
state  operating  condition  to  give  the  linearized  form. 

Hence  it  can  be  written  as 

ax  + a y + d*  a 2 + o 

bx  fcy  d 2 

The  perturbed  variables  are  Ax,  AY,  AZ,  etc.  and  the 

partial  differentials  ( , bz,  etc.)  that  form  the 

« Ty 

coefficients  of  the  perturbed  variables  are  evaluated  at 
steady  state  operating  conditions.  This  technique  was  also 
followed  by  whalley  (2)  in  modeling  the  same  boiler  plant. 
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E.  BOILER  ANALYSIS 

The  boiler  studied  in  this  analysis  is  a Foster  Wheeler 

D-type-  marine  boiler;  which  is  an  oil  fired, two-drum, natural 

circulation  unit  having  a rated  output  of  20,800  lbs/hour  at 
2 

350  lbs/in  gauge,  1200  F.  A cross  sectional  arrangement  of 
the  unit  is  shown  in  Figure  I.The  following  assumptions 
apply  to  the  physical  simplifications  in  each  of  four 
sections. 

1 . Superheater 

(a)  The  inertial  effects  of  superheated  steam  are 
neglected. 

(b)  The  superheater  tubes  are  assumed  to  be  a single 
capacitance  with  restriction  on  the  drum  side  and  another 
restriction  on  the  load  side. 

(cj  Desuperheaters  are  not  considered. 

2.  Downcomer  riser  loop 

/ 

(a)  Only  natural  circulation  exists. 

(b)  No  boiling  takes  place  in  the  downcomers. 

(c)  Vapor  and  liquid  velocities  in  the  riser  are 
identical. 

(d)  Heat-transfer  rates  to  the  boiling  liquid  from  the 
tube  walls  are  proportional  to  the  cuoe  of  temperature 
difference  between  the  wall  and  the  liquid. 

(e)  Steam  quality  is  uniform  in  the  riser. 

(f)  Liquid  temperature  is  always  the  same  as  the 
saturation  temperature  corresponding  to  drum  pressure. 

(g)  Downcomer  liquid  temperature  is  the  same  as  the 
drum  liquid  temperature. 
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3.  Drum 

(a)  There  is  no  temperature  gradient  across  the  drum 
vapor  phase^and  the  temperature  is  always  the  saturation 
temperature  corresponding  to  drum  pressure. 

(b)  The  liguid  phase  has  no  temperature  gradient  other 
than  across  a very  thin  boundary  layer  at  the  drum  surface. 

(c)  Evaporation  or  condensation  rate  in  the  drum  is 
proportional  to  the  difference  of  liguid  and  saturation 
temperatures. 

(d)  Feedwater  temperature  is  assumed  to  be  constant. 

(e)  Liquid-level  changes  due  to  bubble  formation  in  the 
drum  are  neglected. 


4 . Gas  Path 

(a)  The  air -fuel  ratio  is  assumed  to  be  constant. 

(b)  Temperature  of  combustion  gas  entering  superheater 
is  proportional  to  the  firing  rate. 

(c)  Watarwalls  are  lumped  with  the  riser-banks. 

(d)  The  heat-transfer  rate  at  each  tube  bank  is 
determined  by  the  tube  wall  temperature  and  the  average  gas 
temperature. 

(e)  Inertia  of  the  hot  gases  is  neglected,  that  is, 
velocity  changes  take  place  instantaneously. 

(f)  Delays  due  to  the  heat  capacitance  of  the  hot  gases 
are  neglected,  that  is,  temperature  changes  take  place 
instantaneously  in  combustion  gases. 

(g)  All  heat  transfer  is  due  to  turbulent  convection 
and  radiation. 

The  following  steps  are  taken  in  developing  the 
eguations  for  each  of  the  four  sections.  A simple  schematic 
diagram  of  each  sub-section  is  included  at  the  beginning  of 
each  part  of  the  analysis.  A brief  statement  of  the  physical 
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situation  is  included  under  the  headings.  The  dummy 

coefficients  such  as  a,  b , C,  D , were  used  for 

1 1 11  ii 

convenience^ and  their  values  are  declared  in  the  Appendices. 

The  quantities  such  as  5 , h , "h  , etc.  are  the  steady  state 

B S B 

values  of  H , ti  , h , respectively.  These  values  can  be 
3 S B 

evaluated  from  the  original  unperturbed  nonlinear  equations 
by  setting  all  derivative  terms  to  zero  and  solving  for  the 


5 . Superheater  continuity 


Mass  flow  in  - Mass  flow  out  = Bate  of  change  of  enclosed 
oass 

* i(LsAsA) 
at 

Ihe  perturbed  equation  is  : 

iWs-AWs  = LsAsSi/>s 

or, 

M'V  + W a,si/>  (i) 

waere  S is  the  Laplace  variable  denoting  differentiation 
with  respect  to  time 

6 . Superheater  energy  balance 

Heat  input  from  tube  wall  + Heat  in  incoming  steam  - Heat 
with  outgoing  steam  = Rate  of  change  of  internal  energy 

VV,A-WA  - 7^L*V.K*) 

The  perturbed  equation  is  *, 

Since,AhR«o,  A and  L^AgSA/^  -AWR-AWS 
After  substitution  and  rearrangement: 

(LS  - AS.ATs  (*) 

or. 


a+S  ats  + a9AwB  + a6  Aas  - <*7At8 
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7.  Hea t conduction  across  superheater  walls 

Heat  input  to  wall  - Heat  output  from  wall  = Rate  of  change  of 
energy  stored  in  superheater  wall 

Q - GL  * -(\CT  ) 

9S  S dt  S s sw' 

The  perturbed  equation  is  ; 


ag9S  - Afts  - MsCsSAT 


sw 


or. 


VSs*  Vas=atoSATsw 


8.  Heat  input  to  superheater  wall 


(3) 


The  following  empirical  formula  includes  the  cooling  effect 
if  air  supply  differs  from  ideal  air  supply  H (whalley  (2)) 

U.  k nr 


0-6 


The  perturbed  equation  is 


Aa 


UV&Tj  + 2K^2SC*)4tl1 


0.6 


9s 


w 


or. 


aiiAa9s  * ai2Awf  + ai3(^T9S- AtsJ + a4dwa 


C4) 
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9.  Heat  transfer  across  superheater  wall 


For  turbulent  gas  flows: 


Ihe  perturbed  equation  is  : 


• |t.sMtS,-Maw8  + KsWb°'‘(4Tm-4ts) 


* °i,(AT,  - AT, ) 


10.  Heat  transfer  from  combustion  gas 
Average  gas  temperature  at  superheater  wall  is  T where 

' 9s  - TC r— 

^9Wf 

and  C is  related  to  the  heat  capacitance  of  the  combustion 
9 

gas  by  ; 

C3  -c«(l  +P 

where  C is  constant 
c 

Thus  ; 


ags  - 2CcCwf^a)(Tc-T  ) 


The  perturbed  equation  is  ; 


A%&  S 2C^S  + Wa)^Tc  + C\-T90^+(Tc-T,s)AWa-(w^Vla)ATJ 


Bearranging ; 


or, 


-L  iis-  (V i»)ATc+(Te-T„)illf  + (Te-f|J*W<i; 


a A&  - a,a  ATe  + a aw<  + a aw  At 

18  3$  19  c 20  t 2i  a 22u  gs 


1 1 . Superheater  moaentu a 
Neglecting  the  inertia  of  the  steam  ; 

(p>-V\-(l4)A>A  -» 

since  : w6  -/„ASV6 

(p6-  p).ihf*a  )fi  -o 

Substitute  ; l L& 

U *“Vs2 

The  perturbed  eguation  is  ; 


AP_ 


- A Ps  , 2 J AW  - 4 I*  A/3 

7sA2'B 


Hearranging  ; 


2 h AWb  + A Ps  -j&  ^ APb  =,  A PB 
A A 


or. 


Cl  A*.  + a,A  AP3  + AA  “a  APfc 

2 3 B 24  3 28  J B 26  6 


23 


(Wf+Wa)ATgs 

(6) 


(?) 


i 


12.  Riser  continuity  equation 

Mass  flow  into  riser  - Mass  flow  out  = Rate  of  change  of 
enclosed  mass 

\(Lv-V)(/+/,)/2] 


j r 

W - W - — 

" dt 


The  perturbed  equation  is  : 

AW,-AW  - 

since 

1 = 1 + C1'*) 

y°  A A 


then. 


V ’Ak-k)***7T&A 


If  the  starting  point  of  evaporation  in  the  riser  tub.e  is  ^ 
above  the  water  drum  exit  where  ; 

yj  - i - CWw 


then, 

A D - -CAW, 

Ay.  ^ S Aj°  can  be  neglected  since  >j  is  small 

Substituting  Ayo  and  A)j  and  rearranging 


o-sAyL^O  (4-  - ^^ShX  + 0-3  AyLyX  Kft  j 


S A P- 


* 0-5  ^SW)C  SAW, 


= Aw,  - Aw 


or- 


a 7 s A X + ajfts  APRta29SAWw  _ a30Aw*  + a3lAw 


(8) 
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Therefore  : 


After  perturbation  and  substitution  for  AyO  as  done  before 
in  the  Riser  Continuity  Equation,  the  perturbed  equation  is; 


and  is  small  enougn 
to  be  neglected 


h*  * {k  ' J.)  (t;  + f ) -f  L\ 


cc.-i 


- r -* 


Therefore  • 


AP,-  Apft  ~ bvAW  + arSAW  + bKAX-  AWw 

^ Ay-Iw 


or. 


a, 4 p»  + a„A p6  - G3tiw*a„Siwta36ix  + a37Aw 


C9) 


14.  Hear  balance  across  riser  tube  wall 

Heat  from  hot  gas  to  riser  tubes-Heat  from  tubes  to  mixture 
=Hate  of  change  of  internal  energy  stored  in  tubes 

~ a&  ’ dt^  M6CbTRw) 

The  perturbed  equation  is  •, 

- A0S  - M6C,SAt8„ 

oc. 


4 


a ao. 

as  b 


\oS*TBW 


(10) 
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15.  Hi sec  heat  balance  equation 

Heat  input  to  risers-Heat  output  from  risers 
=Sate  of  change  of  internal  energy 

From  the  continuity  aguation  ; 

| [o-sA,(L,-<iHy>+y>1<)]  . w „-w 

After  substitution  and  perturbation  ; 

Ww^kwt>+(^wn-fi)AWw'W  Ah  « 0.5AvLv(/>+yOjs  Ah 

since  ; 

k • kttB+  xkf9 

If  changes  in  latent  heat  Ah^are  neglected , then 

Ah  . A^  + k^&x 

From  steam  tables  ; 

* kcAPb 

A k Ms  a ATW 

Therefore  ; 

0.5  A,Lr(/a^.)KtSaPBt0.SAaL,(/./.)kf,S0X-Aao 

= (k-t-  kJiWj-WK.APj-Sk^ix  + Vi.AT,, 

or. 


«4iS^6 


+ a.,sAx+ail  a a 

42  40  B 


<* A' w + a A P + a A x +a  Atw 

44  * 45  6 46  47  " 


(il) 
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16.  Riser  heat  transfer  equation 

An  eapirical  formula  that  accounts  for  turbulent  gas  flow 
is  adopted  and  gives  t 

V = cWv-tI) 

The  perturoed  equation  is  j 

Wj  * 

♦ aw«  * 4 m v-v-ri  at,,,) 


<xc 


Rearranging 


0.6 


A V *«•**  * V a-VVT&,  + Aw 

"f 

4.  2^<&y(wac-Wft)  Aw  +4Kf3AT 
f " ’a.  + ^ v 'ga  AT9» 


W 


- 4K^Tm  ATbw 


hence  ; 

AQ>"(  '+4^v)AT< 


^0.6 


r 


-(V"|0*  AT0W 


O.fcKga  (T,»  -Taw)  a m. 
+ - °*4-  » 

Vf 

g K»y(wac  - Wa)  Aw 
— 2 
"ac 


or, 


a4a  A ®>< 


+ a49 


a50AT5W  ^a5i^Wf  +aSi4Wa  (12) 
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17.  Riser  heat  transfer  to  boliing  fluid 

The  heat  transfer  rate  is  assumed  to  be  proportional  to  the 
cube  of  temperature  difference  between  wall  and  mixture. 
Continuing  to  follow  Whalley(2): 


- Ka  1 Tb»  ' Tb  ) 


The  perturbed  equation  is  : 


a63  ^ a&  =*  ^54  ATBW  4 ass  ^T6 


C 13 ) 


18.  Downcomer  momentum  equation 

£ Forces  due  to  pressure  difference-frictional  loss 
♦gravitational  head-entrance  loss  = Inertia  force 


Hence  ; 


(P*“P")A«>-  ^CV«)/wAb+-AAo(L*>+2)-  TzP w\ 

U B * V 

d.  ( '-D  fji 


since 


-AAbVrt 


j.  2 

(P B - Pw  ) - 4.0  L - Jil 


W„ 


Vw  2^AyOw  dt  <3A|> 


( p-)-  (^7  + 0-5J 


3A*; 


vt,  4./>  L„  . -ti  i-CW») 
dt 
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The  perturbed  equation  is  ‘ 

ApB  - Ap»  - (b.ta.S^w.-L.A/o. 

where  : 

b«  s 


Since  flow  is  incompressible  ( y>w  is  constant  both  in  time 
and  space), then  _ q 

Therefore  : 

aS6  APa  + ag7  APW  = a58^ww  + asaSAWw  (14.) 

1 9 . Drum  steam  mass  balance  equation 

£ Vapor  rates  to/from  drum  = Rate  of  change  of  vapor  mass 
in  drum 
Thus  , 

Evaporation  rate  ♦ vapor  rate  from  risers-steam 

rate  from  drum  = — ( VBPa) 

dt  J 

Therefore  : 

*e  ♦ XW  - Wft  . ^ Vft) 

The  perturbed  equation  is  : 

AWe  +XAW4  WAX- AWr  - VRS  A/^  4-y5BSAVB 
Substituting  for  ; 

A Va  3 -A  Ay 
A fh  * krAP6 
AWfi  - KeCATw-ATB) 


f 

V 


_2 

ab 


fft  '-D 


+ 0.5 
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gives  ; 


va^BSAP6  -/a ASAy  _ Ke(ATw-AT6)  + x^w  + WAX- AWa 


or. 


a6oSAPB  + a6is  Ay- a62(ATw-ATB)+a63Aw  + a6^AX+a  aw&  (15) 

20.  Drum  liquid  heat  balance  equation 

£ Heat  to  and  from  drum  with  liquid  = Rate  of  change  of  the 
internal  energy  of  the  liquid 
Hence  ; 

Internal  energy  of  riser  liquid  + heat  transfered  in 
feedwater  - heat  taken  in  by  downcomer  liquid  - latent  heat 
released  with  evaporating  liquid  =^(drum  liquid  internal 
energy)  dt 

Therefore  : 

CvU-X)wTB+CiWiTi-CDWwTw-WeKaw  = i(CDMTj 

at 

But  Cy.  , Ci  J Qb  - 1 

The  perturbed  equation  is  ; 

( 1-  X ) T6  A W + (i-x)wat&- wTBAx  + wtATi■^T£Awt.-w*AT* 

-:^'"  AWW-  we  aK8w-Rbw  AWe  - MSATw+ThSAM 

Substituting  for  ; 

4We  » Kc  ( ATW  - ATft) 

A hRy)  . Kc  A Pft 

AM  =,  a/v* a y 

and  I.  #»  constant, so  AT.  * o 
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Thus: 


U-X)fB  AW  + (l-X)WATa-WTBAX  +TiAWi-  WwATw-twAW,-WcKcAP6 

-Krv,Kc(ATw-  ATft)  - M SATW  + A/wTwSAy 
or, 

a66  aw  + a67  atb  + a&s ax  + a6g  aw^  +•  a70ATw  +a71  aww  + a72  apb 
+ ar5(  ATw-  Atb)  - a74SATw  +-a  „SAy  C 16) 

21.  Dr u m liquid  mass  balance  equation 

v Liquid  input/output  rates  = Rate  of  change  enclosed  mass 
of  liquid 
Hence  : 

feedwater  input  rate  ♦ riser  liquid  input  rate-downcomer 
liquid  rate-evaporation  rate  = d ( M ) 

dt 

Therefore  ; 

W;+  ( L-X)  W - Ww  -Wc  s 4(M) 

dt 

The  perturbed  equation  is  ; 

AW;  -I-  (1- X)  AW  - W AX  - AWw-  AWtf  - SAM 

Substituting  for  : 

A We  =*  ^e  (ATW-ATR) 

AM.  Ay°w  Ay 

gives; 

-Ke  ATB  + A/wSAy  - (i-X)  AW-WAX-AWw-KeAT^  + AWt 

0*, 

a« ^ * a„ s 4 y * as,iw+  a,sAX+a.„4vvvvwi  (t7) 
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22.  Equation  from  throttle  valve 

WS  * *1  K,  Ps  >TS) 


The  perturbed  equation  is  : 


C.  STATE  VARIABLE  EQUATIONS 


The  following  procedure  is  an  algebraic  method  to  reduce 
the  twenty-four  governing  equations  to  ten  state  variable 
equations, that  is,  in  standard  matrix  form  of  X(t)  = AX(t)  ♦ 


BU(t) 

Here  X (t)  is  the  matrix  of  system  state  variables  and 
U (t)  is  the  matrix  of  system  inputs. 

From  (1) 


4A 


a.  AWft  +.  a AVI, 
i a 2 & 

A w f l2  Aw 

u3  B a3  s 


or, 


From  (2) 


AA  = biAwB+ 


a+ ATS 


Mws  ~ afeAa  + a At, 


(L.l) 


From  (5) 


AT, 


Aw  _ cu  ag,  + At. 

a*  6 a4  s a4  s 


A a . is?  aws  + “a  at  - at 
6 a is  Qi5  a15 


Ci-2) 


(i.  3) 


Substituting  Q from  (1.3)  in  (1.2)  and  rearranging  ; 


ATa 


[ aU  ai7  ^ 

At  _ 

4 

Ull 

<31 

au 

Ql6 

\Q*‘  ft» 

s 

a4 

a 4. 

cJ 

Q6'  al7 
a4  ,al5 
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ATS  * 


From  (3) 


a 10  s ats.w  * A e?.  +•  a3  A Qs 

At  . = ®lAa  + «i  AQ 
iio  3s  a io  6 


From  (4) 


AV^AWf  + °T  V^4T-+a‘*iw“ 


From  (6) 


A a . = ±12  AT  + AW}  + ii.1  Awa  + at 

°Lia  c au  aia  ai8  3 

Setting  (1.6)  = (1.7)  and  rearranging'. 


4V- 


f<hz_  OacA  ( gl%  • 

Aw,  4-  Aiii-L  ATsa 

/<*13  _ ^aa \ / Qaa  A 

\ttu  ai4j  iau  axaj 

(n  ) f 

Lilliil/  AWsl+  Ul8— - Atc 

/Oi3_aM,\  / £13  _ 

^axi  dxa/  [an  4x1j  J 


(i-4) 


(1.5) 


(i-  6) 


Cl  7) 


( l.  a) 


Substituting  AT  3 Cl  AWr  from  (24)  into  (1.9)  and 

102  1 

rearranging  ; 

/aii  at82  _ ^12  + ®£o  ^ / fl13  \ I a2j  \ 

AT  . . l ~ a“  «»0  4W  + M_  AT  AWa 

I3  /cu*  a21  A * (^13  ^aa \ /aia  a« \ 


bit  uis 


au  axa 
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AT?S  * b7iTSH-b«AWa  (1'9) 

Substituting  AT^5  from  (1.9)  into  (1.6)  and  rearranging  . 

A' V - (; + *£) Aw<  - fc) 4T- + (a“- 

or, 

&a3S=  bs  AW(  + bi0AT„  + baiwa  (1.10) 

Substituting  A a^,  from  (1.10)  and  A&$  from  (1.3)  into 
(1.5)  and  rearranging  ; 


auu  i.  ea*.*.  au  y , 

s"  ^ aw  J ^ [ a t0  alc-aLSj  s*  \ aio  / 


+(V_^±Uv,  - 


au  ' ai-8 


a„  . a. 


a io  • a. 


Ats>  . buAWf  V b13AT3,(.b14AW(x+  b„  aw„  - blt  Ats  (i-lO 


From  (17) 


_ At  +.  ^52  Aw  +-  Ax 
a7,  B a77  a77 


+.  AWw  Atw  + A Wi 

a77  U77  a77 


(M2) 
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Substituting  AT8  from  (19)  into  (1.12) 

A V s _ S6  V AP*  + AW  AX 

' a77  s a77  tt77 

4_  ill  AWw  + °lli  ATW  * ^ AW( 

a77  Q77  a77 


or, 


Ay  = - b7APR  + bu  AW  + bt9  AX+  b20  AWw  + b2lATw  4 b22 Awj 
From  (16) 


ATW 


T.  = - ili.  A }♦  ist  AW  + £*1^  47.  ♦ fu.  « 

w a7+  Q.74.  a 74.  fl7+ 


i»  AW;  * (i2ll i)  A7„  + ilL  AW  I + 

a74  a74  '*■ 


t 

Substituting  Ay  from  (1.13) 
and  rearranging  • 


AT«  , 


^a7S-bi7  a9;-a07  Q97  - ^74  1 A P 4.  ( 

t _ - I ' B In 

a. 74  a 74.  o74 


and  ATb  from  (19)  into 


_ A 75 . bia\ 


a 74 


. f _ a75  • I A X + ( — - 
1 - a 74  / l a74. 


a,7S  bso 


1 74 


74- 


A Wu 


f a7fl  t a73 

a 7+ 


ars  1 bji 


At,+ 


69  «75 


®7o  • bi.  ^ 


•74. 


tt74 


a 


74 


AW: 


or. 


AT,  < b^Ap^.  b244N  + t2sAx  + b.6A«w+  b^  + bj.AW; 


(L- 13) 


(1.14) 

.14) 


(1.15) 
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From  (15) 


ap  _ _ fli1  Ay  atw  _ In  at.  aw  + £m  a*  + i*£Aw  ( 1 . 1 6 ) 
ft  a60  • a.6o  o60  B a60  »6o  o-bo  n 

Substituting  Ay  from  (1.13)  into  (1.16)  and  rearranging  : 

a,.,  (°M  _ aa-bt4) 


APn 


f<hi\r_  M Ap  +(fi3  _ _^L_Ii*Uw  [^4  _ 

^ a60  a60  ) & ^a60  ab0  ) ^ a6« 


C 1 6o 


AX 


AWw  + ( _^L_  fiiibi)  atw 
I J 1 d- to  dto  J 


iii^sUw;  +/  ^]lw 
“U  I ^6.) 


or. 


APa  - b23  A Pft  + b3c  Aw  + bsi  Ax  - b32  Aww 


+ ATw  " ^+Awt  f b36  A wf 


(i  l 7) 


From  (11) 


AX  - - A p . - 513  Aa  + fli  + A P 

a+2  ft  <x42  6 a4a  a42 


f AX  + £il  AT, 
a42  a42 


(1.18) 


From  (13) 


a a 


1*4. 


L S3 


AT 


BN 


“S3 
fl-ff  3 


AT, 


ES 
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Substituting  ATB  froin  (19) 


A 1. 


AT  + °”  a»;  AP 

« 8"  r, _ _ B 


C t • 13 ) 


a53  a53 

I 

Substituting  A Pft  from  (1.17)  and  A Ql&  from  (1.19)  into 
(1.18)  and  rearranging  ; 

A fi  \ ( (\  jl  « li  \ 

AW 


d,A.t.b»a  ®43»®SS*^St^  a n 


Ax  3 I ^l8  - a^-Dg? 

la42  a42 


*»3»aS8,Ct97  ] ^ p _ 

a42'  aB3  J B 


_46 

\a4.1 


a4l-  ^31^  + / «44  + a4l-  b32 


a+2 


a 42 


(X  42 


a42 


AV 


'5te-f«^4T 


‘4i 


0-4-2 


W 


[a 1 Aw- 

f a4b  ■ ^.35  ^ 

Aw 

r a 4.t  • a 34  \ 

a42  j 

a42 

“ 6 ~ 

a42-  ^53  I 

/ 

l i 

^ / 

BW 


or, 

AX  , bjtAPB  - bJ7AW  * b3,AX  + b„  Au, 

+ b40  AT.  f b4l  AW;  - b42  AWft-  b43  AT„.  a 20) 


From  (8) 


A w AX 

d 2g 

AP.  +iAS  AW,  + ^Ai 

AW 

(1  21) 

a29 

a^s 

a19  * aZ9 

0 

Substituting  AX 

from 

• 

(1.20)  and  A PB 

from 

(1.17)  into 

(1.21)  and  rearranging  ; 
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b 

27  .16  , 

^23- b29 j 

a29 

<*29  j 

®2«b3i 

a 23 

u29  j 

®27  • b4fl 

328  • ^33 

a29 

Q-z 3 y 

/ fl27  • b42 

a2«b35\ 

4f,s  + 


ai’  b«  _ a>»- t 2i‘\  4W 

a29  Q.  29  a29 1 


■»'  °32  + U30  _ Q^y  D33 
a 29  &29  ^23 


a29  & 23 


AWw  « - b44  APB  + b4s  AW  - b46  AX  + b47AWw  - b+s  AT* 
+ b49  AWj+  bSo  A Ws  + b5i  ATftw 


(i-22) 


Substituting  AWw  from  (1.22)  into  (14)  and  rearranging  ; 


A Pw  ,(  a5«^a53-  bi7 


l57  } 

\ 


)AW"-(---^7;bt4JAp^(^)4w 


f.  , \ 


-|  °1<l  | A X - ATU  + / *■— — 1 AW. 

l / i J t / 


a-M'  bso  AW  .).  --_39_bl  | AT 


or, 

AW  * b6o  Aww  + b6t  APfl  + b62  Aw  - Ax  - b^  atw 

+ b63AW.  + b66  AW  ± b6  AT  ( i-  2b) 

° f»W 

From  (10) 

at  = Aft,  + Afi  (1  26) 

aW  r a*.  ». 
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Prom  (12) 


a*»  . a50  4T>>  * a„  AW}  * as,  4w4-a4s  4TJo 


(1.2r) 


Since  AT  > AT.S  from  (21)  , then  substituting  AT 

« * 3 ® 

From  (1.9)  into  (1.27)  and  rearranging  j 

Al„  - AT,„  + (.“»>-  “»»  b*  ) AW, 

o <14$  \ a 4.,  / l a 4 8 


f ^43  ^7 

a4S 


AT, 


sw 


( i 28) 


Substituting  A S.gB  from  (1.28)  and  A GLft  from  (1.19)  into 
(1.26)  and  rearranging  ; 


AT 


T - ( Qaa-  a3Q  + q39  • Q-S4 

\ a.4.0 • 0.4.3  04.3.  a33 


1T-+&J(as‘-vb‘)n 


q33  **49  W 


AT 


sw 


or. 


+ f a33  asa-  Q9.7_j  A 

I a4a-aS3  ) 


AThw  = b6s  ATftw  + b69  AWf  4 b70  Awa-  b7i  ATsw  + b72  APb 

Prom  (7) 


(1.29) 


AW  = if!  APfi  _ l£±  A P .^iA/3 

a23  ftfs  5 ai3  y R 


(1.30) 
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Substituting  A Ps  from  (2  3)  and  A/>&  from  (20)  into  (1.30) 


or, 

• bra  APft  - ^74  Ar&  _ b7S  Aj>s  (iai) 

Prom  (18) 

AW c _ AXV  + ^-5  APS  + Jhi.  At  Cl- 32) 

s ~ a93  a93  a93  s 

Substituting  A Ps  from  (23)  into  (1.32) 

Ax,  ,fWi«  ^ £»Ut £ + AA 

l O93  <^93/  4 93 


or, 

^Ws  ■ b76  AXV  + b77  ATs  +.  b78  (i.39) 

.Substituting  AWa  from  (1.31)  and  A W$  from  (1.33)  into 
(1.1)  and  rearranging  ; 

V*  - ( b*  ■ b7,  - k.  ■ • b « ) A/*  + ( ba  ■ b„  - K ■ k„ ) ATS 

+ ( *V  b*>  ) AP,  + C b,  -b,t  ) Axv 

or. 

Vs  - Cn  «A  ♦ Cu  A-s  * Clt  4P„  + Du  AXV  (A) 

Substituting  AWft  from  (1,31)  into  (1.4)  and  rearranging  : 


45 


AT  a b^-b75  A/>s  + ( b3  + b4  - b74,  ) AT3  “ bs  ATSW  - b4  b73  A PB 


or. 


ATS  = ^21  Ay^s  •+  ^22  ATs  + C23  AT3w  + C 23  A F&  ( ) 

Substituting  AW^  from  (1.31)  into  (1.11)  and  rearranging 
ATsw  a - bis,b75  Ay>s  - ( bi5.  b74  + bi6)  ATS  -i-  bi3  ATSW 

+ biS  . b73  A PB  + bu  AWj  + bl4>  AW a 


or. 


Atsw  * ^ 3i  Ay°s  + d32  ATS  + Cam  ATSW  + A PB 

+ &32AWf  + D33^a 


(C) 


Substituting  A WB  from  (1.31)  into  (1.20)  and  rearranging  • 
AX  - b-V2‘b7a  A/s  + b42‘  b74  ATS  + b3S  AX  - b37  AW 
+ b39  AWW  ~ b43  A 1 BW  + ( b3fe“  b4-2  ‘ b7a)  A P& 

i bfa  ATW  + b4l  AWi 


or. 


AX  .CW4/S  tC+jAT,.  + C*4AX  + C45AW  + C46  AW. 
tc47  AT,„+<:4,AP#+C4S  AT*  4- t>44  AW; 


(D) 
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Substituting  A from  (1.31)  into  (1.25)  and  rearranging  ; 


A 


AW  - - ^66-  bra  A/s  - b64  . b7+  Ats  - b63  A X + b62  AW  + b60  AWw 
+ b67  atbw  + ( b66-  b73  + b6i  ) APa  -b6+ATw4  b65  AW( 

or, 

0 

A w s d8l  AyO  + cs2  ATS  + ds+ Ax  + dga  AW  + C56  AW, 

+^57  A76W  +^5®  A Pft  + d5g  ATW  + DS4  AW*  (E) 

Substituting  A WR  from  (1.31)  into  (1.22)  and  rearranging  ; 

* 

AWw  - - bsa  . b75  2y>s  _ b5o-  b74  Ats  - b46  AX  + b4S  A W 
+ b47  AWw  + b51  ATftW+(b50.b73-b44)  AP6 
- b48  atw  + ,b^3  A Wj 
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— - 


The  state  variable  equations  from  (A)  to  (J)  were  rewritten 
in  matrix  notation  as  shown  in  Figure  4.  The  state  equations 
were  soxved  for  various  percent  step  changes  to  input 
variables  ( AXy  , AWj ) , using  the  IBM  simulation  language 
CSMP  - III  as  shown  in  the  next  section. 
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Figure  4 - STATE  VARIABLE  MATRIX 


III.  COMPUTER  RESULTS 


From  Appendix  A , it  is  seen  that  the  a coefficients  are 
in  terms  of  steady  state  values.  These  values  were  • 
determined  from  the  data  in  Whalleys  thesis  (2)  which  are 
repeated  in  Appendix  E.  A Fortran  IV  program  to  find  the 
a,b,C,D  coefficients  which  appear  in  the  CSMP  program  is 
shown  in  Appendix  F,and  the  calculated  values  are  shown  in 
Appendix  G.  Appendix  H is  the  IBM  simulation  language  CSMP 
program  that  was  developed.  Only  input  variables  were 
changed  for  each  computer  run.  For  example,  results  for  a 5% 
step  change  of  A xv , a 10%  step  change  of  AWjand  a 10%  step 
change  of  Awx  were  calculated.  The  output  response  curves 
were  plotted  using  the  NPS  CALCOMP  PLOTTER. 


A.  TRANSIENT  RESPONSES  FOLLOWING  A STEP  CHANGE  IN  THROTTLE 
SETTING  OF  5* 


At  a particular  steady  state  operating  condition  the 
throttle  valve  is  suddenly  opened  while  still  keeping  the 
previous  air  and  fuel  flow  rates  and  feedwater  flow  rate 
constant.  Because  of  bubble  formation  in  the  steam  drum 
liquid, the  steam  drum  water  level  (FIGURE  9)  swells  for 
about  40  seconds  and  then  shows  a steady  falling  in  level  at 
higher  flow  rates  of  steam.  The  steam  drum  pressure  (FIGURE 
8)  immediately  begins  to  fall.  Following  these  sudden 
increases,  the  steam  mass-flow  rates  from  the  steam  drum  and 
the  superheater  (FIGURES  10  and  11)  show  slight  declines 
because  of  the  effect  of  the  decrease  in  steam  drum 
pressure.  The  superheater  outlet  pressure  (FIGURE  10)  acts 
in  a similar  fashion  to  the  steam  drum  pressure,  that  is  ,it 
shows  a steady  decline  to  a new  steady  state  condition  but 
it  is  slightly  lower  than  the  steam  drum  pressure  because  of 
system  pressure  drops.  The  superheated  steam  temperature  and 
superheater  wall  temperature  (FIGURES  5 and  6)  both  decline, 
alt-aough  the  latter  temperature  is  slightly  lower.  Tae 
riser  tube  wall  temperature  (FIGURE  8)  and  steam  drum  liquid 
temperature  (FIGURE  9)  also  decrease  due  tc  the  throttle 
change.  The  superheated  steam  density  (FIGURE  5)  shows  a 
rising  characteristic  to  counteract  the  loss  in  superheater 
pressure  and  temperature  and  the  gain  in  steam  mass-flow 
rate  from  tne  steam  drum  (FIGURE  11).  The  riser  and 
downcomer  flow  rates  (FIGURE  7)  are  reduced  to  a lower  value 
but  steam  quality  (FIGURE  6)  is  increased  to  a higher  value 
than  previously  to  preserve  the  energy  balance. 
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Figure 


5 - SUPERHEATED  STEAM  DENSITY  < A/>$  ) VS  TIME  AND 

SUPERHEATED  STEAM  TEMPERATURE  < AT.  > VS  TIME 


5.  CO 


Figure  6 - SUPERHEATER  WALL  TEMPERATURE  < AT„  > vs  TIME 
AND  QUALITY  OF  MIXTURE  LEAVING  RISER  < AX  > VS  TIME 


54 


5-00 


-d.od  b. oo  J2.00  ie.no  24.03  30.00 

TIKE  rxio  ] 1 


Figure  8 - KISER  TUBE  RALL  TEMPERATURE  < ) vs  TIME  AMD 

STEAM  DRUM  PRESSURE  < A PB  > VS  TIME 
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Figure  9 - STEAM  DRUM  LIQUID  TEMPERATURE  < Alw  > VS  TIME 
AND  STEAM  DRUM  HATER  LEVEL  < Ay  ) VS  TIME 
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Figure  11  • STEAM  MASS-FLOW  HATE  FROM  STEAM  DRUM  TO 
SUPERHEATER  < > VS  TIME 
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B.  TRANSIENT  RESPONSES  FOLLOWING  A STEP  CHANGE  OF  1056  IN 
FUEL  FLOW  RATE. 


The  effects  due  to  a change  in  the  fuel  flow  rate  are 
less  than  those  due  to  a change  in  the  throttle  valve 
setting.  However, the  system  requires  a longer  period  of  time 
to  reach  a new  steady  state  condition.  Because  the 
combustion  rate  is  increased,  the  evaporation  rate  is  also 
increased  making  both  the  steam  mass-flow  rate  at  the 
superheater  (FIGURE  17)  and  the  steam  mass-flow  rate  from 
the  steam  drum  (FIGURE  18)  increase  in  a similar  manner. 
Initially,  density  in  the  risers  falls  quickly  in  response 
to  the  increased  firing  rate  which  causes  a rise  in 
superheated  steam  density  (FIGURE  12)  . An  increase  in  steam 
flow  from  the  steam  drum  without  an  increase  in  the 
feedwater  flow  rate  results  in  a drop  in  steam  drum  water 
level  (FIGURE  16). The  steam  drum  pressure  (FIGURE  15)  and 
superheater  outlet  pressure  (FIGURE  17)  rise  monotonically 
as  time  increases.  Also  superheated  steam  temperature 
(FIGURE  12)  , superheater  wall  temperature  (FIGURE  13)  , 

riser  tube  wall  temperature  (FIGURE  15)  and  steam  drum 
liquid  temperature  (FIGURE  16)  all  increase  with  increasing 
fuel  flow  rate.  The  abrupt  fall  in  riser  mixture  flow  rate 
(FIGURE  14)  is  probably  due  to  the  sudden  drop  in  the  level 
of  the  initial  density  change  in  the  riser. 
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Figure  12  - SUPERHEATED  STEAM  DENSITY  < A/>s  ) VS  TIME  AND 
SUPERHEATED  STEAM  TEMPERATURE  < AT-  ) VS  TIME 


Figure  13  - SUPERHEATER  HALL  TEMPERATURE  < ^Ts»(  > VS  TIME 

AND  QUALITY  OF  MIXTURE  LEAVING  RISER  < AX  » VS  TIME 
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Figure  14  - RISER  MIXTURE  FLOW  RATE  < AW  ) yS  TiBE  AND 
DOWNCOMER  LIQUID  FLOW  RATE  < AW*  > VS  TIME 
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AND  STEAM  DRUM  WATER  LEVEL  < Ay  > VS  TIME 
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C.  TRANSIENT  RESPONSES  FOLLOWING  A STEP  CHANGE  OF  10%  IN 
FEEDWATER  FLOW-RATE 


The  main  effect  of  increasing  feedwater  flow  rate  while 
maintaining  throttle  setting  and  air-fuel  flow  rate  constant 
is  an  increase  in  the  steam  drum  water  level  (FIGURE23). 
Because  feedwater  is  added  into  the  steam  drum,  steam  drum 
pressure  (FIGURE22)  suffers  a small  overshoot  followed  by  a 
decline.  As  a result,  a drop  in  both  steam  mass-flow  rate 
from  the  steam  drum  (FIGDRE25)  and  steam  mass-flow  rate  at 
the  superheater  outlet  (FIGURE24)  occurs. Flow  around  the 
riser-downcomer  loop  (FI3URE21)  is  only  sligntly  affected. 
Since  the  firing  rate  is  unalter ed, the  riser  tube  wall 
temperature  (FIGURE  22)  and  the  steam  drum  liquid 
temperature  (FIGURE23)  decline.  The  superheated  steam 
temperature  (FIGURE19)  and  the  superheater  wall  temperature 
(FIGURE20)  rise  because  of  a drop  in  steam  mass-flow  rate 
from  the  steam  drum. The  quality  of  steam  (FIGURE20)  shows  a 
general  lowering.  The  superheated  steam  density  (FIG  UR  El  9) 
shows  a decline  whicn  is  similar  in  form  to  the  drop  in  the 
superheater  outlet  pressure. 
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Figure  20  - SUPERHEATER  WALL  TEMPERATURE  < AT$W  ) VS  TIME 
AMD  QUALITY  OF  MIXTURE  LEAVING  RISER  < AX  ) VS  TIME 
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Figure  22  - RISER  TUBE  HALL  TEMPERATURE  < AT#W  ) vs  TIaE 
AND  STEAM  DRUM  PRESSURE  C A P&  ) yS  TIME 
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Pigure  24  - SUPERHEATER  OUTLET  PRESSURE  < APS  > VS  TIME 
AND  STEAM  PLOW  RATE  AT  SUPERHEATER  OUTLET  < > VS  TIME 
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IV.  CONCLUSIONS 


This  analysis  is  only  the  initial  step  toward  a better 
understanding  of  a naval  steam  generating  plant  both  from 
the  dynamic  and  from  the  control  points  of  view.  The 
analysis  concerns  a boiler  system  which  is  divided  into  four 
parts.  The  major  difficulty  in  the  analysis  is  the  fact 
that  the  whole  system  is  very  complex  and  contains  numerous 
variables  which  are  extremely  unwieldy  to  manipulate.  Large 
numbers  of  relationships  are  non-linear  so  it  is  necessary 
to  generate  a linear  form  by  the  methods  of  perturbation 
theory. After  setting  up  the  state  variable  equations,  the 
IBM  simulation  language  CSMP-HI  was  used  to  solve  tor  the 
open  loop  transient  response  characteristics  of  the  state 
variables  due  to  variations  in  input  variables.  Throughout 
these  calcu lations, all  tne  variables  appearing  in  the 
transient  response  curves  are  not  absolute  values. 
Bather, they  represent  a small  incremental  change  from  some 
particular  steady  state  operating  condition.  Also, it  is  to 
be  noted  that  trends  in  some  of  these  results  differ 
markedly  in  form  from  those  appearing  in  Hhalley  (2) . Future 
studies  will  be  undertaken  to  determine  the  causes  of  these 
differences.  Meanwhile, the  benefit  one  may  derive  from  thi's 
basic  analysis  is  its  use  in  the  design  and  construction  of 
a new  type  of  boiler  controller,  that  is,  a multivariable 
control  system. 


76 


APPENDIX  A 


LIST  OF  COEFFICIENTS  A 


a = 1 


a = - 1 


S ~ VsMs 


a = L A 
4 S S 


a = h -h 
s SB 


a = - 1 


a = -W  C 
7 S ps 


a = 1 


a = - 1 


a = M C 
io  s S 


a = 1 
i i 


77 


_ 0-6  f \ 

a!2 3s  TswJ 

_0.6 

a = Kgj 
13  3 7 


L = 2 K Cflac- WaJ 
14  w1 

" oc 


a = 1 


1 3 


0.8  Ks  fT  'l 

= ■_ V 'aw  - 


16  w 


& 

o.  b 


1 7 


K H 
S B 


« 2Cc 


a = (H  *W  ) 
i»  fa 


20 


21 


22 


(T  -T  ) 
C gS 


(I  -T  ) 
C gS 


-<W  *v  ) 

f a 


a = 2f  JSft. 

23  s A2 


a = 1 

24 


2 S 


7» 


t _ri£. 

TS  *2 


78 


n c 

B B 


a41  = 0‘5  Ay  Ly  [p  + y°w)  Kc 


a42  = 0.5A.L,C/5 


a = - 1 

4 3 


a = <h  -h) 

♦♦  HD 


a = -HK 
♦5  C 


a = - Hh 

♦ ft  fg 


a = u 


S»Wf  -4K,t38 


» =-s**r-4MB3w 


0-6  *9f>(T3»-T6>)) 
- 0,4 
Wf 


2 ( *<ac  ~ W g ) 


a =1 

S3 


84 


b 

12 

_ 

a io 

• 

b 

1 3 

_ /a,,  bto 
\ a ie 

b 

14 

Qa.  bn 

a io 

’ 

b 

15 

®9  • a16 

Q-ii-  al5 

b 

16 

aB  • ai7 

aw  ais 

b 

1 7 

_ a7fe-a97 

<177 

h 

II 

Is 

\ 

u 

14 

a 77 

b 

QiS3 

M 

1 9 

a 77 

b 

20 

a9o 

a 77 

L 

b 

2 1 

_ tt91 
®77 

h 

0-92 

XJ 

22 

*77 

Qa -^tr 

•^lfl  ■ aiS 


a63  ^46 


a59- V 
O-sj 


^59.^49 


fl39  ' ^ 30 


aS9-*>51 

a3r 


Q32- ^32 

v 

a3« 

a3S 

Vl  - 

a3I.  ^33 

a35 

a36 

°32 -^34 

_ £at 

Q3J '^58  + a36 

0-1  wr  tt.H 


Vie 


a3g*  ^37 


_ a3£-k3 


90 


APPENDIX  D 
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b2  7=(  A70  4-W?  )/47  ♦-<  A 7r>*rin  /A7  ♦ ) 
&2  3*(AoO/a7^)-(  a?5*6  2 2/A  •'♦) 


B30=l  A63/A6J)-  ( A6l*HL'1/A&3  ) 
b 2 1 = ( A0A/A6  J)  - ( A 6 l*P  lr/A63) 


635=  A i>  5 / AoJ 

A 45 /A  V>1-  ( AU*tJ2r-/A4‘>  ) -<  A A3*  AO  5*A9'W4'V  A *3  ) 


r3S=(  AA0/A-. : ) - ( AAl*b  U/ A4'>  ) 
h2  0*(  »w'*/ aas if(  4i i*;: ’ VA<r  1 
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R43=(A+3*\54>/ 

ls44=(  AZ  7*6  5 j/A2  3)  + ( A?"*B25/A?^  ) 


6 A 0=1  AZ7-s.-i3  -)/M’  7 )Hi'2a^.M/A',c) 

t»A/'  = ( A?  3*i  3 2 /A  2 ■:  J +(  A30/A29  J-(  A ?7*93  9/425  ) 


t;Ac  = < A2»«H34/«2,7)-(  A?7*3A1  7A?3  ) 
i * •-»  4 ’ / A iV?fi*H“.S/A  2 


o5l=A2  7*rt-»3/A2l> 

BS»  2- ( A53/A57)  + ( A 5 <J  *l»  A’’/ A5~  ) 


nj4  = Aj?*Ljib/A:  7 

h ‘5  5=  m 5 9* c4t//  uc’  •’ 


h o J=  ( A 2 £ rt  > 5 Z / A ?.  f> ) - ( A ? 7 / 4 7 3 ) 

i / >.  ■<  *>  y i * » i • - ■>  . . i c * j k ^ r i 


66  2=  (A3  2 *6  3 4/  A 2 *i ) - t A 3 A / A?  •>  ) 

I?  6 5=(  A i <*  3 6 6/  a J 3 i + ( A 3 £ / A*  ‘J  ) 


( 4 1 = ‘ R 7 ? 
Ch2«342-*h  7 A 
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3cX«5hl3;.  =* 


6 -t  X • b h L t 
7SX, 


= .r-l  = .?,5X,5HC:3  BtfL5.e//t 

— c 1 c c c v c l r -)  1 _ i~  i r-  n # / 


.^X,5i-:3rf  s,ff  l5.°f5X,5wr^l  = ♦ F 1 5 • 9//  t 
. b*  , 5t-  r.V'+  »,E1  5.8f  bX,5hr.45  =,FI5.9//, 


~ tC  '4,r  QS  t 

IC  -5X  tc  3C  ,u  v,C9  4,C  >5  ,C16  9 8 . C< S , C 10 4 #C4 05 , C l'J6  , C IC6  ,C.IC9 

roVIATllt-  » ti  X » b r*  f.  t ‘j  = ,M5.3,5Xt5HC*6  = » El  f>  .9 , 5X  » 5hC67  *,F15.9//, 


»•  - ’ » ' >'  V . T V V • V I f > ■ \ J P *.  • C f ^ f I ’/  r „ " f u A t u r / f 

2SX  1 5riC  7 7 «,Elb.9«  5X,  5HC7?  * , r 1 5 . P , 5 X , 5HC  9 1 =,E15.?//» 
3 VX»  5'«.3Z  »»«:lS.t>,3X,3>-C3A  =*clE»a»5X*5  wf  9 9 *,ei5.e//. 


•t  V A • • I)  - t l A v • " * -J  A » 3-v.  5 " * t r I.  7 • ^ t f 3«U6V  =«C  / t 

5C  X * 5.-*C  <:•  ^ =,E13.3,5X,‘7FC9k  « , f 1 E .b  , 5 X , 5HC.  ?A  = ,E15  •£//♦ 
A c a . 9m  < ,•  = . 1 1 s . i i y. . Ai-r  t t s.pm.e.sjy.ijumts.tK.a//. 
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APPENDIX  G 

FORTRAN  IV  PROGRAM  OUTPUT 


a i o")s  j .:o)Q)::i  3:> 
*73  = 0 .1415274 Cc  ->32 
WTT'  = 0 • ♦ 2 -J  J 'J  0 0 G L JO 


Cll  = -C.  l.243t0,Z0r:  J2 

;r.  34 


C4?  =-  0 . (j 0 7 j 2 2 l C - 0? 


6 


1 i ci-  i~. 353: conob  05 
t>74  =C.  13333340F  0 ) 

' > 7 ? =""' . 1 5 'Tf-  — 37“ 

ci'  =-c.n4Ut,5or-ci 


0.13  19 4440 F 03 

T.rcC^o^ir^r 

0.1315921 06-03 
0. 5 1. c 563406  01 


C44  =-0 . 4'76ft2"2  26  01 


= C.27Ccjr,3uCE-0r> 


CL' 

C5  7 = 

C 6 


=-G  • t>6  7 35  5 006 -0 1 

C 7 £ 

5 = 0. 255? 11 106-33 

m 

= O.Bi3333-'Cfc  03 


■0  • 3 005 1 5 c Ofc-0  i 
■0. 375  04ft  IC2-.J3 


C fc  4 = 0.  44097750P  OS 


CSS  = -0.?’Cr>820CF  00 


Cl3ft=-0.72Cl25l0F-03 


032  = 0 . 3 4 C 9 4 ? 5 OF  01 
u54  =-0.5:5154202  02 


C. 14577520c  C4 


3. 9*2562706  05 


= 0.3298 61 3 OF  06 


C 3 5 = C.ll 41CS70E  03 


c 


cio4=-o.i2eeei5CE  co 


Cl 0 3=- 0.3  45 9 4270 E- 06 


032  = 0.6576360 06 -01 
064  =-C. 774236606  01 


r 
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APPENDIX  H 


CSMP-HT  PEOGfi  AM 


"!R  HO SP  =C  1 il '*0  T S ♦ C l:*rP^  U1 

t'PHns*!  NTiiftL  I ' C * , )R,-Cr.  C) 


OT  S = '•  NT.J.KL  ( " C 2 t U T S 0 J 

kISWD=C2  l.*D^H.'S  +C:  ?.f  IT  S *i".  33  * rTc.  w K i - * - ow  ; 5 *.  c ^ r-  - 


V I j A = 1 i y I yKL  \ ^ • I V*  ’I 

j x u».c  *♦  1 * 0 k k i s * c a ti  v • j r .$  f c 4 -♦ 

i.r  /.o*i\nw  *.r  • t . A «•#  -a  . • 


nx  = iNTviKt_<  rc, 't.rxi:) 

rwD»C5l*OkHlSfC5«’*i>T!>*C:>'»*''X+r55**'»  . .. 


l'r  1 ri,n  1*11  1 »r  n t •••.  n 

l)lnW®  IfiT  ijKl.  ( f C t tl'V,  <j) 


UTf}w=!NTGkl  ( ’ C 7 » n f oV-,:j ) 
nP  ilO  = CB  l *CRHO$  + C*‘  l*H  TS  ♦r.p  '<*r  > «-f.  ' * • >V  +r  a <•  *!)•»  • , ff  C>  o*r'r>.»  f (•  ]cj  ; T»  ... 


tHIS  PAQS 

raoMCOPx 


is  BtST  9UU*111  FRAfillOW* 


PAGE  HE’C.HT=3»U'nTh  = S 
Pa  G L XlTPLUT 


LA 6 cl.  OT?  Vi  T'  VF 
OUTPUT  T : M E • 0 T S U 

i f act  •>  t c we  -r  * ■ 


OUTPUT  THE. OX 
LABEL  OX  VS  T M t 


i r v/  i i 

LABEL  Ow  VS  T vc 
OUTPUT  TVVE.  D»V! 


k a d hi.  u w 'a  Vi  i ' v t 
OUTPUT  T'SF.OTBV. 

I A J El  ,"\T  Al.  \ I C T * r 


OUTPUT  THE.GPB 
L AC  C.L  UPB  VS  T v,  c 


output  thf.dv^s 
l A Ctrl  J‘*s  Vi  T'.  Mr 


hm  r u i » 41  it 

LABEL  OWB  VS  t ; F 

t mep  f :f.t- ^ioo.u.r-i.noL-1  = 2.0 .p’iof’  = -.c 
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